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I. INTRODUCTION 

In the last few years a great deal of attention has been paid to the study of interacting 
quantum fields out of equilibrium. There are, in fact, many interesting physical situations in 
which the standard S-matrix approach cannot give sensible information about the behavior 
of the system, because it evolves through a series of highly excited states (i.e., states of finite 
energy density). 

As an example consider any model of cosmological inflation: it is not possible to extract 
precise predictions on physical observables without including in the treatment the quantum 
back-reaction of the field on the space-time geometry and on itself |]T]-§[. 

On the side of particle physics, the ultra-relativistic heavy-ion collisions, scheduled in the 
forthcoming years at CERN-SPS, BNL-RHIC and CERN-LHC, are supposed to produce 
hadron matter at very high densities and temperatures; in such a regime the usual approach 
based on particle scattering cannot be considered a good interpretative tool at all. To 
extract sensible information from the theory new computational schemes are necessary, that 
go beyond the simple Feynmann diagram expansion. The use of resummation schemes, 
like the Hartree-Fock (HF) [^Q approximation and the large N limit (LN) , or the Hard 
Thermal Loop resummation for systems at finite temperature (HTL) [0, can be considered a 
first step in this direction. They, in fact, enforce a sum over an infinite subset of Feynmann 
diagrams that are dominant in a given region of the parameter space, where the simple 
truncation of the usual perturbative series at finite order cannot give sensible answers. 

Quite recently HF and LN have been used in order to clarify some dynamical aspects 
of the large A^ 0^ theory, reaching the conclusion that the non-perturbative and non-linear 
evolution of the system might eventually produce the onset of a form of non-equilibrium 
Bose-Einstein condensation of the long-wavelength Goldstone bosons usually present in the 
broken symmetry phase p|- p!0|] of the model. Another very interesting result in ^ con- 
cerns the dynamical Maxwell construction, which reproduces the fiat region of the effective 
potential in case of broken symmetry as asymptotic fixed points of the background evolution. 

In this article we present a detailed study, in finite volume, of dynamical evolution out 
of equilibrium for the $^ scalar field. More precisely, we determine how such dynamics 
scales with the size of the periodic box containing the system in the case of uniform back- 
grounds. This is necessary to address questions like out-of-equilibrium symmetry breaking 
and dynamical Bose-Einstein condensation. 

We apply two of the non-perturbative methods mentioned above, namely the Hartree- 
Fock approximation and the large A^ expansion. 

In section |l| we define the model in finite volume, giving all the relevant notations and 
definitions. We also stress the convexity of the effective potential as an exact result, valid 
for the full renormalized theory in any volume. 

In section |T^ we derive the large A^ approximation of the 0(A^)— invariant version of 
X{(f)'^y model, according to the general rules of ref. |]ll|. In this derivation it appears evident 



the essential property of the N ^ oo limit of being a particular type of classical limit, so 
that it leads to a classical phase space, a classical hamiltonian with associated Hamilton's 
equations of motion [see eqs. (|3.15|) , (|3.16|) and ( p. 171) ]. We then minimize the hamiltonian 



function(al) and determine the conditions when massless Goldstone bosons (i.e. transverse 
fluctuations of the field) to form a Bose-Einstein condensate, delocalizing the vacuum field 



expectation value (cfr. also ref. |T2|). This necessarily requires that the width of the zero- 



mode fluctuations becomes macroscopically large, that is of the order of the volume. Only 
when the background takes one of the extremal values proper of symmetry breaking the 
width of the zero-mode fluctuations is of order L^/^, as typical of a free massless spectrum. 
The study of the lowest energy states of the model is needed for comparison with the 
results of the numerical simulations, which show that the zero-mode width o"o stays micro- 
scopic (that is such that ctq /volume— > when the volume diverges) whenever it starts from 
initial conditions in which it is microscopic. Our results, in fact, show clearly the presence of 
a time scale r^, proportional to the linear size L of the system, at which flnite volume effects 
start to manifest. The most remarkable consequence of the presence of such a scale is that 
it prevents the zero mode amplitude to grow macroscopically large. This result contradicts 
the interpretation of the linear late-time growth of the zero-mode width as a Bose-Einstein 
condensation of Goldstone bosons |3-[T0[ . 



On the other hand we do flnd that the size of the low-lying widths at time r^ is of order 
L, to be compared to the equilibrium situation where they would be of order L° in the 
massive case or of order L^/^ in the massless case. Perhaps the denomination "microscopic" 
should be reserved to this two possibilities. Therefore, since our initial condition are indeed 
microscopic in this restricted sense, we do observe in the out-of-equilibrium evolution a 
rapid transition to a different regime intermediate between the microscopic one and the 
macroscopic one characteristic of Bose-Einstein condensation. 

At any rate, when one considers microscopic initial conditions for the choice of bare mass 
which corresponds to broken symmetry, the role itself of symmetry breaking is not very 
clear in the large N description of the out-of-equilibrium dynamics, making equally obscure 
the issues concerning the so-called quantum phase ordering [^. This is because the limit 
A^ ^ oo is completely saturated by gaussian states, which might signal the onset of symmetry 
breaking only developing macroscopically large fluctuations. Since such fluctuations do not 
appear to be there, the meaning itself of symmetry breaking, as something happening as 
times goes on and accompanied by some kind of phase ordering, is quite unclear. Of course, 
in this respect the main limitation of our approach, as well as of those of the references 
mentioned above, is in the assumption of a uniform background. Nonetheless, phenomena 
like the asymptotic vanishing of the effective mass and the dynamical Maxwell construction, 
taking place in this contest of a uniform background and large N expansion, are certainly 
very signiflcant manifestations of symmetry breaking and in particular of the Goldstone 
theorem which applies when a continuous symmetry is broken. 



To gather more information on these matters, we consider in section |^ a more elab- 
orate type of time-dependent Hartree-Fock (tdHF) approximation, which generalizes the 
standard gaussian self-consistent approach to non-gaussian wave-functionals; in fact, one 
might envisage the possibility that, while gaussian fluctuations never become of the size of 
the volume, non-gaussian fluctuations do grow in time to a macroscopic size. We derive 
therefore the mean-fleld coupled time-dependent Schroedinger equations for the modes of 
the scalar fleld, under the assumption of a uniform condensate, see eqs ( ^4.5| ), ( ^^ ) and ((4.7]). 



A signiflcant difference with respect to previous tdHF approaches [0] concerns the renor- 
malization of ultraviolet divergences. In fact, by means of a single proper substitution of 
the bare coupling constant Ab with the renormalized one A in the Hartree-Fock hamiltonian, 
we obtain completely cut-off independent equations (apart from the corrections in inverse 



power of the cutoff which are there due to the Landau pole). The substitution is introduced 
by hand, but is justified by simple diagrammatic considerations. 

One advantage of not restricting a priori the self-consistent HF approximation to gaussian 
wave-functionals, is in the possibility of a better description of the structure of the vacuum 
in case of broken symmetry. In fact we can show quite explicitly that, in any finite volume, 
the ground state the zero-mode of the field is concentrated around the two vacua of the 
broken symmetry, driving the probability distribution for any sufficiently wide smearing of 
the field into a two peaks shape. This is indeed what one would intuitively expect in case 
of symmetry breaking. On the other hand none of this appears in a dynamical evolution 
that starts from a distribution localized around a single value of the field in the spinodal 
region, confirming what already seen in the large A^ approach. More precisely, within a 
further controlled gaussian approximation of our tdHF approach, one observe that initially 
microscopic quantum fiuctuations never becomes macroscopic, suggesting that also non- 
gaussian fiuctuations cannot reach macroscopic sizes. As a simple confirmation of this fact, 
consider the completely symmetric initial conditions (0) = (0) = for the background: in 
this case we find that the dynamical equations for initially gaussian field fiuctuations are 
identical to those of large A^ [apart for a rescaling of the coupling constant by a factor of 
three; compare eqs. ( |3.15| ) and ( |4.28| )], so that we observe the same asymptotic vanishing of 



the effective mass. However, this time no interpretation in terms of Goldstone theorem is 
possible, since the broken symmetry is discrete; rather, if the width of the zero-mode were 
allowed to evolve into a macroscopic size, then the effective mass would tend to a positive 
value, since the mass in case of discrete symmetry breaking is indeed larger than zero. 

On the other hand, also in the gaussian HF approach, we do find that the size of the 
low-lying widths at time tl is of order L. We then discuss why this undermine the self- 
consistency of the gaussian approximation, imposing the need of further study, both analyt- 
ical and numerical. 

Finally, in section |V| we summarize the results presented in this article and we sketch 
some interesting open problems that we plan to study in forthcoming works. 

II. CUTOFF FIELD THEORY 

We consider the scalar field operator in a /}— dimensional periodic box of size L and 
write its Fourier expansion as customary 



X 



L-^''Y.<P>^e 



ik-x i\ 

' Vk 



with the wavevectors k naturally quantized: k = {27c/L)n, n G Z-^. The canonically conju- 
gated momentum tt has a similar expansion 

k 

with the commutation rules [0^ , i^-k'] = i^kJ ■ The introduction of a finite volume should be 
regarded as a regularization of the infrared properties of the model, which allows to "count" 
the different field modes and is needed especially in the case of broken symmetry. 



To regularize also the ultraviolet behavior, we restrict the sums over wavevectors to the 
points lying within the D— dimensional sphere of radius A, that is /c^ < A^, with TV = KL/2tv 
some large integer. Clearly we have reduced the original field-theoretical problem to a 
quantum-mechanical framework with finitely many (of order M^~^) degrees of freedom. 

The d)^ Hamiltonian reads 



H=- I d^x 
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where m^ and Ab should depend on the UV cutoff A in such a way to guarantee a finite 
limit A — ^ oo for all observable quantities. As is known [^jQ) this implies triviality (that 



is vanishing of renormalized vertex functions with more than two external lines) for D > 3 
and very likely also for D = 3. In the latter case triviality is manifest in the one-loop 
approximation and in large— iV limit due to the Landau pole. For this reason we shall keep 
A finite and regard the 0^ model as an effective low-energy theory (here low-energy means 
practically all energies below Planck's scale, due to the large value of the Landau pole for 
renormalized coupling constants of order one or less). 

We shall work in the wavefunction representation where ((y9| \1/) = ^['(v?) and 

d 

(0o^)(</') = ^^-^W) , K^)(¥') = -^J^'^{v) 

while for A; > (in lexicographic sense) 

1 1 f d d \ 

(0±fc^)(<^) = ^i^k±t ^-k) ^(y^) , (vr±fc^)(v?) = -^ I -^^ ± d—J ^(<^) 

Notice that by construction the variables (fk are all real. Of course, when either one of 
the cutoffs are removed, the wave function "^{f) acquires infinitely many arguments and 
becomes what is usually called a wave functional. 

In practice, the problem of studying the dynamics of the 0^ field out of equilibrium 
consists now in trying to solve the time-dependent Schroedinger equation given an initial 
wavefunction '^{'^^t = 0) that describes a state of the field far away from the vacuum. 
By this we mean a non-stationary state that, in the infinite volume limit L ^ oo, would 
lay outside the particle Fock space constructed upon the vacuum. This approach could be 
generalized in a straightforward way to mixtures described by density matrices, as done, for 



instance, in |16-18|. Here we shall restrict to pure states, for sake of simplicity and because 
all relevant aspects of the problem are already present in this case. 

It is by now well known ||T^ that perturbation theory is not suitable for the purpose 
stated above. Due to parametric resonances and/or spinodal instabilities there are modes 
of the field that grow exponentially in time until they produce non-perturbative effects for 
any coupling constant, no matter how small. On the other hand, only few, by now standard, 
approximate non-perturbative schemes are available for the 0^ theory, and to these we have 
to resort after all. We shall consider here the time- dependent Hartree-Fock (tdHF) approach 
(an improved version with respect to what is presented, for instance, in 0]) and the large N 



expansion to leading order. In fact these two methods are very closely related, as shown for 
instance in [|l^], where several techniques to derive reasonable dynamical evolution equations 
for non-equilibrium (p'^ are compared. However, before passing to approximations, we would 
like to stress that the following rigorous result can be immediately established in this model 
with both UV and IR cutoffs. 



A. A rigorous result: the effective potential is convex 

This is a well known fact in statistical mechanics, being directly related to stability 
requirements. It would therefore hold also for the field theory in the Euclidean functional 
formulation. In our quantum-mechanical context we may proceed as follow. Suppose the 
field is coupled to a uniform external source J. Then the ground state energy Eq{J) is 
a concave function of J, as can be inferred from the negativity of the second order term 
in AJ of perturbation around any chosen value of J. Moreover, Eq{J) is analytic in a 
finite neighborhood of J = 0, since J(f) is a perturbation "small" compared to the quadratic 
and quartic terms of the Hamiltonian. As a consequence, this effective potential KflF(<^) = 
Eq{J) — Jcj), (j) = Eq{J) = {(f))Q, that is the Legendre transform of Eq{J), is a. convex analj^ic 
function in a finite neighborhood of = 0. In the infrared limit L -^ oo, Eq{J) might 
develop a singularity in J = and Ves{4>) might flatten around = 0. Of course this 
possibility would apply in case of spontaneous symmetry breaking, that is for a double-well 
classical potential. This is a subtle and important point that will play a crucial role later on, 
even if the effective potential is relevant for the static properties of the model rather than 
the dynamical evolution out of equilibrium that interests us here. In fact such evolution is 
governed by the CTP effective action |n||20| and one might expect that, although non-local 
in time, it asymptotically reduces to a multiple of the effective potential for trajectories 
of (f){t) with a fixed point at infinite time. In such case there should exist a one-to-one 
correspondence between fixed points and minima of the effective potential. This is one of 
the topics addressed in this paper. 

III. LARGE N EXPANSION AT LEADING ORDER 

A. Definitions 

In this section we consider a standard non-perturbative approach to the 0^ model which is 



applicable also out of equilibrium, namely the large N method as presented in [^. However 
we shall follow a different derivation which makes the gaussian nature of the N —>■ oo limit 
more explicit. 

In the large N method one generalizes the 0^ model by promoting the single real scalar 
field to a A^— component vector cf) of scalar fields, in such a way to ensure 0{N) symmetry. 
This corresponds to the hamiltonian 



2 



4L^ 

k ki,k2,k3,k4 



where the space integration and the sum over wavevectors are hmited by the infrared and 
ultraviolet cutoff, respectively, according to the general framework presented in section |l|. 

It is known that this theory is well behaved for large A^, provided that the quartic coupling 
constant A is rescaled with 1/A^. For example, it is possible to define a perturbation theory, 
based on the small expansion parameter 1/A^, in the framework of which one can compute 
any quantity at any chosen order in 1/A^. From the diagrammatic point of view, this 
procedure corresponds to a resummation of the usual perturbative series that automatically 
collects all the graphs of a given order in 1/A^ together p. Moreover, it has been established 
since the early 80's that the leading order approximation (that is the strict limit A^ — > cxd) is 



actually a classical limit |TT[, in the sense that there exists a classical system (i.e., a classical 
phase space, a Poisson bracket and a classical Hamiltonian) whose dynamics controls the 
evolution of all fundamental quantum observables, such as field correlation functions, in the 
N ^ oo limit. For instance, from the absolute minimum of the classical Hamiltonian one 
reads the energy of the ground state, while the spectrum is given by the frequencies of small 
oscillations about this minimum, etc. etc.. We are here interested in finding an efficient and 
rapid way to compute the quantum evolution equations for some observables in the N ^ oo 
limit, and we will see that this task is easily accomplished just by deriving the canonical 
Hamilton equations from the large N classical Hamiltonian. 



Following Yaffe [11|, we write the quantum mechanical hamiltonian as 

H = Nh{A,C) (3.1) 

in terms of the square matrices A, C with operator entries {vuk is the canonical momentum 
conjugated to the real mode ipk) 

A ^ ^1 

^kk' — ~j^'^k' ^k' 1 ^kk> — l^'^k' '^k' 

These are example of "classical" operators, whose two-point correlation functions factorize 
in the N ^ oo limit. This can be shown by considering the coherent states 



^2,q,p(9^) = C{z) exp 



'Y.Pk-V'k-^Y.^kk'i^k-Qk)- i^fik'-Qk' 

k kk' 



(3.2) 



where the complex symmetric matrix z has a positive definite real part while p^ and g^ are 
real and coincide, respectively, with the coherent state expectation values of tu^ and (^^. As 
these parameters take all their possible values, the coherent states form an overcomplete set 
in the cutoff Hilbert space of the model. The crucial property which ensures factorization 



is that they become all orthogonal in the N ^ oo limit. Moreover one can show [|ri[| that 
the coherent states parameters form a classical phase space with Poisson brackets 



{lk,Pk'} 



P.B. 



^kk'S'''' , {Wkk' y'Vqq'}p,B. ~ ^kq^k'q' + ^kq'^k'q 



where w and v reparametrize z as z = \w~^ + iv. It is understood that the dimensionality 
of the vectors qk and pk is arbitrary but finite [that is, only a finite number, say n, of pairs 
{(pi , wl) may take a nonvanishing expectation value as A^ -^ oo]. 



Once applied to the classical operators Akk' and Ckk' the large A^ factorization allow to 
obtain the classical hamiltonian by simply replacing A and C in eq. (|3.1|) by the coherent 
expectation values 

{Akk') = Qk ■ Qk' + Wkk' , {Ckk') = Pk ■ Pk' + {vwv)kk' + j{w'^)kk' 

In our situation, having assumed a uniform background expectation value for 0, we have 
Qk — Pk — ^ for all k ^ 0; moreover, translation invariance implies that w and v are diagonal 
matrices, so that we may set 

_ 2 X — ^^ X 
Wkk' — <^k "kk' 5 Vkk' — Okk' 

0"fc 

in term of the canonical couples {ak, Sk) which satisfy {ak , Sk'}p.B. = ^kk'- Notice that the 
(Tfc are just the widths (rescaled by N~^^'^) of the 0{N) symmetric and translation invariant 
gaussian coherent states. 

Thus we find the classical hamiltonian 



hci = ^ipl + ml ql) + -Y, 



sl + {e + ml)al + 



4ai 



^ A. 



4L^ 



{^0 + T.^l)' 



where by Hamilton's equations of motion Pq = ^q and Sk = cFk- The corresponding conserved 
energy density £ = L~^hc\ may be written 



2 2L^ ^ 

(3.3) 

2 
k 



v-^i:(*^^.^ + ^)+W^-.5:). ^-^^^i 



where (f) = L'^^'^Qq and V is the 0(A^)— invariant quartic potential regarded as a function 
of 0^, that is V{z) = ^mlz + |Ab-2^. It is worth noticing that eq. ( |373|) would apply as is to 
generic V{z). 

B. Static properties 

Let us consider first the static aspects embodied in the effective potential Ves{4'), that 
is the minimum of the potential energy V at fixed (j>. We first define in a precise way the 
unbroken symmetry phase, in this large A^ context, as the case when V^fj(0) has a unique 
minimum at = in the limit of infinite volume. Minimizing V w.r.t. ak yields 

4 = , ^ , M^ = m^ + 2 V'{S^ + S) 

"" 2Vk^ + M2 ' ^ ^^ ^ 

2 , \ J,2 , Ab -r-^ 1 

that is the widths characteristic of a free theory with self-consistent mass M fixed by the gap 
equation. By the assumption of unbroken symmetry, when = and at infinite volume M 
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coincides with the equihbrium mass m of the theory, that may be regarded as independent 
scale parameter. Since in the hmit L ^ oo sums are replaced by integrals 



d'^k 



(Ti 



we obtain the standard bare mass parameterization 

ml = m'-\Jn{m\A), Id{z,A)^[ ^f f„ ^ i (3.5) 



fc2<A2 (27r)« 2VWT 



and the renormalized gap equation 



M' = m' + X(f)' + X[lD{M\A)-lD{m\A)\ . (3.6) 

which implies, when D = 3, 

X^ = x(l-^\og^] (3.7) 

with a suitable choice of the finite part. No coupling constant renormalization occurs instead 
when D = 1. The renormalized gap equation ( |3.6|) may also be written quite concisely 





M- 


' "^ , ;i,^ 


A(M) ~ A(m) ' ^ 


in terms of the one-loop running couplings constant 


A(/i) = A 


1 ^ , /" 


-1 

, A(m) = A 



(3.8) 



It is the Landau pole in A(2Ae^^/^) that actually forbids the limit A -^ oo. Hence we must 
keep the cutoff finite and smaller than Apoie, so that the theory does retain a slight inverse- 
power dependence on it. At any rate, there exists a very wide window where this dependence 
is indeed very weak for couplings of order one or less, since Apoie = |mexp(l/2+87r^/A) ^ m. 
Moreover, we see from eq. (p.8|) that for \/A|0| much smaller than the Landau pole there 
are two solutions for M, one "physical", always larger than m and of the same order of 
m + vA|0|, and one "unphysical" , close to the Landau pole. 

One can now easily verify that the effective potential has indeed a unique minimum in 
= 0, as required. In fact, if we assign arbitrary 0— dependent values to the widths ak, 
(minus) the effective force reads 

^^(*.{..W})^M=^- + i:f$ (3.9, 

and reduces to M^ 0* when the widths are extremal as in eq. (|3.4|) ; but M^ is positive for 
unbroken symmetry and so = is the unique minimum. 

We define the symmetry as broken whenever the infinite volume Ves has more than one 
minimum. Of course, as long as L is finite, Ves has a unique minimum in = 0, because of 
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the unicity of the ground state in Quantum Mechanics, as aheady discussed in section P A| . 
Let us therefore proceed more formally and take the limit L — > oo directly on the potential 
energy V. It reads 




V=7TL^.„77TZTn(A:^«^,^ + — 1+^0 +S), S= /,„,.„ tTTZwJ ^. 



fc2<A2 (2tt) 



where we write for convenience the tree-level potential V in the positive definite form V{z) = 
jXh{z + m^/Ab)^. V is now the sum of two positive definite terms. Suppose there exists a 

configuration such that V{(p + S) = and the first term in V is at its minimum. Then this 
is certainly the absolute minimum of V. This configuration indeed exists at infinite volume 
when D = 3: 

0^' = ^, ^' = v\ ml = -X^[v' + h{0,A)] (3.10) 

where the nonnegative v should be regarded as an independent parameter fixing the scale 
of the symmetry breaking. It replaces the mass parameter m of the unbroken symmetry 
case: now the theory is massless in accordance with Goldstone theorem. On the contrary, 
ii D = 1 this configuration is not allowed due to the infrared divergences caused by the 
massless nature of the width spectrum. This is just the standard manifestation of Mermin- 
Wagner-Coleman theorem that forbids continuous symmetry breaking in a two-dimensional 



space-time 22 



At finite volumes we cannot minimize the first term in V since this requires uq to diverge, 

— 2 

making it impossible to keep V{4> + S) = 0. In fact we know that the unicity of the ground 
state with finitely many degrees of freedom implies the minimization equations ( ^.4[ ) to hold 
always true with a M^ strictly positive. Therefore, broken symmetry should manifest itself 
as the situation in which the equilibrium value of M^ is a positive definite function of L 
which vanishes in the L ^ oo limit. 

We can confirm this qualitative conclusion as follows. We assume that the bare mass has 
the form given in eq. ( p. 101 ) and that (p = v'^ too. Minimizing the potential energy leads 



always to the massive spectrum, eq. ( |3.4|) , with the gap equation 

M2 1 1^1 A' 



+ oT3z2 /m-im s^2 (2-^^) 



vr^ 



Ab 2L3M 2L3 ^0 y/k^ + M2 8 

If M^ > does not vanish too fast for large volumes, or stays even finite, then the sum on 
the modes has a behavior similar to the corresponding infinite volume integral: there is a 
quadratic divergence that cancels the infinite volume contribution, and a logarithmic one 
that renormalizes the bare coupling. The direct computation of the integral would produce 
a term containing the M^ log(A/M). This can be split into M^[log{A/v) — log(M/f)] by 
using V as mass scale. The first term renormalizes the coupling correctly, while the second 
one vanishes if M^ vanishes in the infinite volume limit. 
When L — > cxD, the asymptotic solution of ( |3.11| ) reads 



/x\l/3 

M= - L-^ + h.o.t. 
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that indeed vanishes in the hniit. Note also that the exponent is consistent with the as- 
sumption made above that M vanishes slowly enough to approximate the sum over A; 7^ 
with an integral with the same M. 

Let us now consider a state whose field expectation value (f) is different from v"^. If 

— 9 

<p > V, the minimization equations ( |3.4| ) leads to a positive squared mass spectrum for the 
fiuctuations, with M^ given self-consistently by the gap equation. On the contrary, as soon 
as < f^, one immediately see that a positive M^ cannot solve the gap equation 

if we insist on the requirement that aQ not be macroscopic. In fact, the r.h.s. of the previous 
equation is negative, no matter which positive value for the effective mass we choose, at 
least for L large enough. But nothing prevent us to consider a static configuration for which 
the amplitude of the zero mode is macroscopically large (i.e. it rescales with the volume 
L^). Actually, if we choose 

we obtain the same equation as we did before and the same value for the potential, that is 
the minimum, in the limit L —>■ 00. Note that at this level the effective mass M needs not 
to have the same behavior in the L ^ 00 limit, but it is free of rescaling with a different 
power of L. We can be even more precise: we isolate the part of the potential that refers to 
the zero mode width ctq (S' does not contain the Uo contribution) 

1 2 \ 4 1 

if o , /-2 „A1 0"n Ab CTg 



mb 



+ Ab (0' + S' 



+ ^T^ + 



L3 4 L6 8LVg 



and we minimize it, keeping (f) fixed. The minimum is attained at t = <Jq/L^ solution of 
the cubic equation 

Abt^ + aAbt^ - jL-^ = 



-,2 



where a = (j) — f + S' — /a (0, A). Note that Aba depends on L and it has a finite limit in 
infinite volume: X{4> — v'^). The solution of the cubic equation is 

Abt = Xy,{v' - 0') + llL'iv' - 0')]~2 + h.o.t. 

from which the effective mass can be identified as proportional to L~^. The stability equa- 
tions for all the other modes can now be solved by a massive spectrum, in a much similar 
way as before. 

Since ctq is now macroscopically large, the infinite volume limit of the o"fc distribution 
(that gives a measure of the transverse fiuctuations in the 0{N) model) develop a 5— like 
singularity, signalling a Bose condensation of the Goldstone bosons: 

4 = {v' - ^') 6^^\k) + ^ (3.12) 

11 



At the same time it is evident that the minimal potential energy is the same as when cp 
that is the effective potential ffattens, in accord with the Maxwell construction. 

Eq. ( p.l2| ) corresponds in configuration space to the 2— point correlation function 



lim 



(0(^) • 0(y)) 

N 



0^ 






2 ik-{x-y) 



C{4>') 



+ ^D{x-y) 



where A£)(x — y) is the massless free-field equal-time correlator, while 
C(0') = v^ Q{v^ - 0') + 0' 6(0' - v^) = max(t;2, 0^) 



(3.13) 



(3.14) 



This expression can be extended to unbroken symmetry by setting in that case C{4> ) = 4> ■ 

Quite evidently, when eq. ( |3.14| ) holds, symmetry breaking can be inferred from the limit 

Ix — 1/1 — >• cxD, if clustering is assumed p3| , ^ , since ^o^x — y) vanishes for large separations. 

Of course this contradicts the infinite volume limit of the finite-volume definition, = 



lim7v-»oo A^ ^^'^{4'{x)), except at the extremal points <p = v^ . 



-,2 



In fact the L ^ oo limit of the finite volume states with 4> < t; do violate clustering, 
because they are linear superpositions of vectors belonging to superselected sectors and 
therefore they are indistinguishable from statistical mixtures. We shall return in more 



detail on this aspects in section |I^, where a generalized HF approximation is considered for 
A^ = 1. For the moment we can give the following intuitive picture for large N. Consider any 
one of the superselected sectors based on a physical vacuum with = t>^. By condensing 
a macroscopic number of transverse Goldstone bosons at zero-momentum, one can build 
coherent states with rotated (f). By incoherently averaging over such rotated states one 
obtains new states with field expectation values shorter than v by any prefixed amount. In 
the large A^ approximation this averaging is necessarily uniform and is forced upon us by 
the residual 0{N — 1) symmetry. 



C. Out of equilibrium dynamics 

We now turn to the dynamics out of equilibrium in this large A^ context. It is governed 
by the equations of motion derived from the total energy density S in eq. (|3.3[), that is 



C^'0 , ^9 -, CP<^k /, 9 , ^9x 1 



4af 



where the generally time-dependent effective squared mass M^ is given by 

M'^ = m^ + Ab [0' + S - loirn^, A) 
in case of unbroken symmetry and 



M2 = Ab 
for broken symmetry in D = 3. 



-,2 



-t;^ + S-/3(0,A) 



(3.15) 



(3.16) 



(3.17) 
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At time zero, the specific choice of initial conditions for at that give the smallest energy 
contribution, that is 



c^fc = , Of. 



2VW+W 



(3.18) 



turns eq. ( p.l6|) into the usual gap equation ( ^^ . For any value of this equation has one 
solution smoothly connected to the value M = m at (p = 0. Of course other initial conditions 
are possible. The only requirement is that the corresponding energy must differ from that 



of the ground state by an ultraviolet finite amount, as it occurs for the choice (|3.18|) . In fact 
this is guaranteed by the gap equation itself, as evident from eq. ( p.9| ): when the widths a^ 
are extremal the effective force is finite, and therefore so are all potential energy differences. 

This simple argument needs a refinement in two respects. 

Firstly, in case of symmetry breaking the formal energy minimization w.r.t. ak leads 
always to eqs. ( p.l^ ), but these are acceptable initial conditions only if the gap equation 
that follows from eq. ( p.l7|) in the L — > oo limit, namely 



M' = Xy 



4>' 



+ /^(M2,A)-/^(0,A) 



(3.19) 



admits a nonnegative, physical solution for M^. 

Secondly, ultraviolet finiteness only requires that the sum over k in eq. ( |3.9| ) be finite 
and this follows if eq. (|3.18| ) holds at least for k large enough, solving the issue raised in the 
first point: negative M^ are allowed by imposing a new form of gap equation 



M^ = Xy 



4>'~v' + 



1 



1 



E ^1 + td E 



1 



fc2<|Af2 



L^ 



fc2>|J\/2| 



2JA;2-|M2| 



/d(0,A) 



(3.20) 



where all cTk with k"^ < \M'^\ are kept free (but all by hypothesis microscopic) initial condi- 
tions. Of course there is no energy minimization in this case. To determine when this new 
form is required, we observe that, neglecting the inverse-power corrections in the UV cutoff 
we may write eq. (|3.19D in the following form 



A(M) 



<p — V 



(3.21) 



-,2 



•2 and 
M'^ = 0, only provided 0' > v'^. So, in the large A^ limit, as soon as we start with (p" < v"^, 



There exists a positive solution M smoothly connected to the ground state, (f) = v" 



we cannot satisfy the gap equation with a positive value of M^. The situation will be quite 
different in the case of A^ = 1 in the HF approximation (cfr. section P^VB| ). 

Once a definite choice of initial conditions is made, the system of differential equations 
( p.l5| ), ( |3.16| ) or ( p.lTp can be solved numerically with standard integration algorithms. This 
has been already done by several authors [pJPiHi working directly in infinite volume, with 
the following general results. In the case of unbroken symmetry it has been established that 
the (Tfc corresponding to wavevectors k in the so-called forbidden bands with parametric 
resonances grow exponentially in time until their growth is shut off by the back-reaction. 
For broken symmetry it is the region in fc— space with the spinodal instabilities caused by 
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an initially negative M^, whose widths grow exponentially before the back-reaction shutoff. 
After the shutoff time the effective mass tends to a positive constant for unbroken symmetry 
and to zero for broken symmetry (in D=3), so that the only width with a chance to keep 
growing indefinitely is ctq for broken symmetry. 

Of course, in all these approaches the integration over modes in the back-reaction S 
cannot be done exactly and is always replaced by a discrete sum of a certain type, depending 
on the details of the algorithms. Hence there exists always an effective infrared cutoff, albeit 
too small to be detectable in the numerical outputs. A possible troublesome aspect of this 
is the proper identification of the zero-mode width a^. Even if a (rather arbitrary) choice 
of discretization is made where a ctq appears, it is not really possible to determine whether 
during the exponential growth or after such width becomes of the order of the volume. Our 
aim is just to answer this question and therefore we perform our numerical evolution in finite 
volumes of several growing sizes. Remanding to the appendix for the details of our method, 
we summarize our results in the next subsection. 



D. Numerical results 

After a careful study in D = 3 of the scaling behavior of the dynamics with respect to 
different values of L, the linear size of the system, we reached the following conclusion: there 
exist a L— dependent time, that we denote by tl, that splits the evolution in two parts; for 
t ^ tl, the behavior of the system does not differ appreciably from its counterpart at infinite 
volume, while finite volume effects abruptly alter the evolution as soon as t exceeds tl] in 
particular 

• Tl is proportional to the linear size of the box L and so it rescales as the cubic root of 
the volume. 

• Tl does not depend on the value of the quartic coupling constant A, at least in a first 
approximation. 

The figures show the behavior of the width of the zero mode ctq (see Fig. |l]), of the squared 
effective mass M^ (see Fig. ^ ) and of the back-reaction S (see Fig. ^, in the more 
interesting case of broken symmetry. The initial conditions are chosen in several different 
ways (see the appendix for details), but correspond to a negative M^ at early times with 
the initial widths all microscopic, that is at most of order L^/^. This is particularly relevant 
for the zero-mode width o"o, which is instead macroscopic in the lowest energy state when 
(f) < v^, as discussed above. As for the background, the figures are relative to the simplest 
case (j) = = 4>, but we have considered also initial conditions with > 0, reproducing 
the "dynamical Maxwell construction" observed in ref. ^. At any rate, for the purposes 
of this work, above all it is important to observe that, due to the quantum back-reaction, 
M^ rapidly becomes positive, within the so-called spinodal time []8|j9|JT^, and then, for times 
before tl, the weakly dissipative regime takes place where M^ oscillates around zero with 
amplitude decreasing as t~^ and a frequency fixed by the largest spinodal wavevector, in 
complete agreement with the infinite-volume results |Q. Correspondingly, after exponential 
grow until the spinodal time, the width of the zero-mode grows on average linearly with time, 
reaching a maximum for t ^ r^. Precisely, ctq performs small amplitude oscillations with the 
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same frequency of M^ around a linear function of the form A + Bt, where A,B^ A~^/^ (see 
Fig. ^), confirming what aheady found in refs. [§],^; then quite suddenly it turns down and 
enters long irregular Poincare-like cycles. Since the spinodal oscillation frequency does not 
depend appreciably on L, the curves of ctq at different values of L are practically identical 
for t < tl. After a certain number of complete oscillations, a number that scales linearly 
with L, a small change in the behavior of M^ (see Fig. ^) determines an inversion in ctq (see 
Fig. ^, evidently because of a phase crossover between the two oscillation patterns. Shortly 
after tl dissipation practically stops as the oscillations of M^ stop decreasing in amplitude 
and become more and more irregular, reflecting the same irregularity in the evolution of the 
widths. 

The main consequence of this numerical scenario is that the linear growth of the zero- 
mode width at infinite volume cannot be consistently interpreted as a form of Bose-Einstein 
Condensation (BEC) ||^. If a macroscopic condensation were really there, the zero mode 
would develop a S function in infinite volume, that would be announced by a width of the 
zero mode growing to values 0(L^/^) at any given size L. Now, while it is surely true that 
when we push L to infinity, also the time tl tends to infinity, allowing the zero mode to 
grow indefinitely, it is also true that, at any fixed though arbitrarily large volume, the zero 
mode never reaches a width 0(L^/^), just because t^ oc L. In other words, if we start from 
initial conditions where ctq is microscopic, then it never becomes macroscopic later on. 

On the other hand, looking at the behavior of the mode functions of momenta k = 
{2tt/ L)n for n fixed but for different values of L, one realizes that they obey a scaling 
similar to that observed for the zero-mode: they oscillate in time with an amplitude and a 
period that are 0{L) (see fig. |^ and H). Thus, each mode shows a behavior that is exactly 
half a way between a macroscopic amplitude [i.e. 0(L'^/^)] and a usual microscopic one [i.e. 
at most 0(L^/^)]. This means that the spectrum of the quantum fluctuations at times of 
the order of the diverging volume can be interpreted as a massless spectrum of interacting 
Goldstone modes, because their power spectrum develops in the limit a l/A;^ singularity, 
rather than the 1/k pole typical of free massless modes. As a consequence the equal-time 
field correlation function [see eq. ( |3.13| )] will fall off as \x — y\~^ for large separations smaller 
only than the diverging elapsed time. This is in accord with what found in [Q, where the 
same conclusion where reached after a study of the correlation function for the scalar field 
in infinite volume. 

The fact that each mode never becomes macroscopic, if it started microscopic, might be 
regarded as a manifestation of unitarity in the large A^ approximation: an initial gaussian 
state with only microscopic widths satisfies clustering and clustering cannot be spoiled by 
a unitary time evolution. As a consequence, in the infinite-volume late-time dynamics, the 
zero-mode width o"o does not play any special role and only the behavior of cr^ as fc -^ is 
relevant. As already stated above, it turns out from our numerics as well as from refs. 



that this behavior is of a novel type characteristic of the out-of-equilibrium dynamics, with 
(Tfc OC 1/k. 

IV. TIME-DEPENDENT HARTREE-FOCK 

The main limitation of the large A^ approximation, as far as the evolution of the widths 
ak is concerned, is in its intrinsic gaussian nature. In fact, one might envisage a scenario 

15 



in which, while gaussian fluctuations stay microscopic, non-gaussian fluctuations grow in 
time to a macroscopic size. In this section we therefore consider a time-dependent HF 
approximation capable in principle of describing the dynamics of non-gaussian fluctuation 
of a single scalar fleld with 0^ interaction. 

As anticipated in section |, we examine in this work only states in which the scalar fleld 
has a uniform, albeit possibly time-dependent expectation value. In a tdHF approach we 
may then start from a wavefuction of the factorized form (which would be exact for free 
flelds) 



k>0 



(4.1) 



The dependence of ipk on its two arguments cannot be assumed to factorize in general 
since space translations act as 50(2) rotations on ipk and ip-k (hence in case of translation 
invariance ipk depends only on (fil'^f'^k)- The approximation consists in assuming this form 
as valid at all times and imposing the stationarity condition on the action 



Sldt{tdt-H)=0, {■) = {^it)\-\^{t)) 



(4.2) 



with respect to variations of the functions tpk- To enforce a uniform expectation value of </> 
we should add a Lagrange multiplier term linear in the single modes expectations {(fk) for 
k ^ 0. The multiplier is then flxed at the end to obtain {(pk) = for all k ^ 0. Actually one 
may verify that this is equivalent to the simpler approach in which (ipk) is set to vanish for all 
k y^ before any variation. Then the only non trivial expectation value in the Hamiltonian, 
namely that of the quartic term, assumes the form 



d^'x {<P{xY) 



1 

+ 



2\2 



{vl)-Hvt) 



+ 



2LD 



EM + ^-.)')-2((^D + (^^))' 



A;>0 



L^ 



U^l) 



(4.3) 



Notice that the terms in the flrst row would cancel completely out for gaussian wavefunctions 
ipk with zero mean value. The last term, where the sum extends to all wavevectors k, 
corresponds instead to the standard mean fleld replacement (0^) -^ 3(0^)^. The total 
energy of our trial state now reads 



E={H) = lj:{:^ + ik' + ml) A + ^Jd^x {m') 



(4.4) 



and from the variational principle (|4.2|) we obtain a set of simple Schroedinger equations 



idtipk = Hki^k 



(4.5) 



Ho 
Hk 



15' 1 2 2 

■2 9^ + 2"°^° + 4L^ 



Ab 4 



M 



+ 



&" 



d^l 



2^kifl 



^\) 



3Ab 

8L^ 



{^i+^\y 



(4.6) 
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which are coupled in a mean-field way only through 

ul = e + ml + ?>K^,, S, = -^ 5: {^D (4.7) 

-^ g2<A2 

and define the HF time evolution for the theory. By construction this evolution conserves 
the total energy E of eq. (|4.4|). 

It should be stressed that in this particular tdHF approximation, beside the mean-field 
back-reaction term E^ of all other modes on uj\, we keep also the contribution of the diagonal 
scattering through the diagonal quartic terms. In fact this is why S^ has no contribution 
from the fc— mode itself: in a gaussian approximation for the trial wavefunctions ipk the 
Hamiltonians H^ would turn out to be harmonic, the quartic terms being absent in favor of 
a complete back-reaction 

^ = ^k + j^ = JD l^M (4-8) 

k 

Of course the quartic self-interaction of the modes as well as the difference between E 
and Efc are suppressed by a volume effect and could be neglected in the infrared limit, 
provided all wavefunctions tpk stays concentrated on mode amplitudes ipk of order smaller 
than L^/^. This is the typical situation when all modes remain microscopic and the volume 
in the denominators is compensated only through the summation over a number of modes 
proportional to the volume itself, so that in the limit L ^ oo sums are replaced by integrals 

Indeed we shall apply this picture to all modes with k ^ 0, while we do expect exceptions 
for the zero-mode wavefunction ipQ. 

The treatment of ultraviolet divergences requires particular care, since the HF approx- 
imation typically messes things up (see, for instance, remarks in |^). Following the same 
login of the large N approximation, we could take as renormalization condition the require- 
ment that the frequencies ul are independent of A, assuming that m^ and Ab are functions 
of A itself and of renormalized A— independent parameters m'^ and A such that 

^l = e + m' + 3\ [S,],^.,^ (4.9) 

where by [.jfinite we mean the (scheme-dependent) finite part of some possibly ultraviolet 
divergent quantity. Unfortunately this would not be enough to make the spectrum of energy 
differences cutoff-independent, because of the bare coupling constant Ab in front of the 
quartic terms in H^ and the difference between E and E^ [such problem does not exist 
in large A^ because that is a purely gaussian approximation]. Again this would not be a 
problem whenever these terms become negligible as L — i> oo. At any rate, to be ready to 
handle the cases when this is not actually true and to define an ultraviolet-finite model also 
at finite volume, we shall by hand modify eq. ([4.3|) as follows: 



Ab / d^'x {<P{xf) =AL-^ {^D - ^^lf + I E. 



fc>0 



((^i + V^'-.n-2((v^i) + (^^_,) 



+ 3AbL^E2 
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We keep the bare coupling constant in front of the term containing E^ because that part 
of the hamiltonian is properly renormalized by means of the usual cactus resummation ^ 
which corresponds to the standard HF approximation. On the other hand, within the same 
approximation, it is not possible to renormalize the part in curly brackets of the equation 
above, because of the factorized form ( |4.1| ) that we have assumed for the wavefunction of the 
system. In fact, the 4— legs vertices in the curly brackets are diagonal in momentum space; 
when we contract two or more vertices of this type, no sum over internal loop momenta 
is produced, so that all higher order perturbation terms are suppressed by volume effects. 
However, we know that in the complete theory, the wavefunction is not factorized and loops 
contain all values of momentum. This suggests that, in order to get a finite hamiltonian, 
we need to introduce in the definition of our model some extra resummation of Feynmann 
diagrams, that is not automatically contained in this self-consistent HF approach. The only 
choice consistent with the cactus resummation performed in the two-point function by the 
HF scheme is the resummation of the 1-loop fish diagram in the four-point function. This 
amounts to the change from Ab to A and it is enough to guarantee the ultraviolet finiteness 
of the hamiltonian through the redefinition 



Hn — ^ Hn 



A-A 



b 4 



Hk 



Hk + 



3(A - Ab) 



{^l + <^' 



At the same time the frequencies are now related to the widths (v5?./c) by 



UJi- 



M^ 



k' + M'-3XL'^i{^t) + {^U)), k>0 



c^o + 3AL-^(^^) = mt + 3AbS 



(4.10) 



(4.11) 



Apart for 0{L~^) corrections, M plays the role of time-dependent mass for modes with 
/c 7^ 0, in the harmonic approximation. 

In this new setup the conserved energy reads 



E 



fc>0 



[XhL S 



+ |AL 



-D 



fc>0 ^ ^ 



(4.12) 



Since the gap-like equations ( [4.11|) are state-dependent, we have to perform the renor- 
malization first for some reference quantum state, that is for some specific collection of 
wavefunctions ipk] as soon as ml and Ab are determined as functions A, ultraviolet finiteness 
will hold for the entire class of states with the same ultraviolet properties of the reference 
state. Then an obvious consistency check for our HF approximation is that this class is 
closed under time evolution. 

Rather than a single state, we choose as reference the family of gaussian states 
parametrized by the uniform expectation value {(f){x)) = L^^^'^{ipo) = 4> (recall that we 
have {(fk) = when A; 7^ by assumption) and such the HF energy E is as small as possible 
for fixed 0. Then, apart from a translation by L^/'^cf) on y^o, these gaussian ifjk are ground 
state eigenfunctions of the harmonic Hamiltonians obtained from Hk by dropping the quar- 
tic terms. Because of the fc^ in the frequencies we expect these gaussian states to dominate 
in the ultraviolet limit also at finite volume (as discussed above they should dominate in the 
infinite- volume limit for any fc 7^ 0). Moreover, since now 

{vl) = L^'^' + ^, {Vlk) = l^. MO (4.13) 



2uJn 



2uk 



the relation ( [4.11|) between frequencies and widths turn into the single gap equation 



"'--'^M^'^^M'7J?Tw) <"^' 



fixing the self-consistent value of M as a function of 0. It should be stressed that ( |4.11| ) 



turns through eq. (|4.13|) into the gap equation only because of the requirement of energy 



minimization. Generic ipk, regarded as initial conditions for the Schroedinger equations 
([4.5|), are in principle not subject to any gap equation. 

The treatment now follows closely that in the large N section, the only difference being 
in the value of the coupling, now three times larger. In fact, in case of 0{N) symmetry, the 
quantum fluctuations over a given background {(f>{x)) = decompose for each k into one 
longitudinal mode, parallel to (j>, and A^ — 1 transverse modes orthogonal to it; by boson 
combinatorics the longitudinal mode couples to (j) with strength 3Ab/iV and decouple in the 
N ^ oo limit, while the transverse modes couple to with strength (A^ — l)Ab/A^ -^ Xh] 
when A^ = 1 only the longitudinal mode is there. 

As L —>■ cxD, a;| — > /c^ + M^ and M is exactly the physical mass gap. Hence it must 
be A— independent. At flnite L we cannot use this request to determine ml and Ab, since, 
unlike M, they cannot depend on the size L. At inflnite volume we obtain 



M2=m2+3Ab[0' + /D(M2,A)], Id{z,A)= f 

Jk 



d^k 1 



fe2<A2 (27r)^2v/FT^ 



(4.15) 



When = this equation flxes the bare mass to be 

ml = m'^- 3Ab/D(m2, A) (4.16) 

where m = M{(j) = 0) may be identifled with the equilibrium physical mass of the scalar 
particles of the inflnite-volume Fock space without symmetry breaking (see below). Now, 
as in large A^, the coupling constant renormalization follows from the equalities 



M^ = m^ + 3Ab[0' + Id{M\ A) - /^(m^. A)] 
= m^ + 3A 0^ + 3A [/d(M^ A) - /^(m^ A) 



(4.17) 

finite 



A 3A , 2A ,, , 

- = 1-— log^ 4.18 



and reads when D = 3 



that is the standard result of the one-loop renormalization group p3|. When D = 1, that is 
a 1 + 1— dimensional quantum fleld theory, /^(M^, A) — l£){m?^h) is already flnite and the 
dimensionfull couplings constant is not renormalized, Ab = A. 

Again, the Landau pole in Ab prevents the actual UV limit A — i> oo. Nonetheless, 
neglecting all inverse powers of the UV cutoff when Z) = 3, it is possible to rewrite the gap 



equation (|4.17| ) as in eq. (|3.8|) : 



A(M) \{m] 
19 



+ 3 02 (4.19) 



in terms of the one-loop running couplings constant 



A(^) = A 



3A /i 

log — 



Stt 



m 



It is quite clear that the HF states for which the renormalization just defined is sufficient 
are all those that are gaussian-dominated in the ultraviolet, so that we have [compare to 
eq. iH)] 



i^lk) 



1 



2uk 



er^A\ A 



oo 



(4.20) 



If this property holds at a certain time, then it should hold at all times, since the 
Schroedinger equations ( ^.5| ) are indeed dominated by the quadratic term for large Uk and 



UJt 



k"^ + const + 0{k ^) as evident from eq. 



Thus this class of states is indeed closed 



under time evolution and the parameterizations ( [4.16| ) and ( |4.18| ) make our tdHF approx- 
imation ultraviolet finite. Notice that the requirement ([4.201) effectively always imposes a 
gap equation similar to eq. ( |4.14D in the deep ultraviolet. 

Another simple check of the self-consistency of our approach, including the change in 
selected places from Ab to A, as discussed above, follows from the energy calculation for 
the gaussian states with {(p{x)) = (j) introduced above. Using eq. ( [4. 4]) and the standard 
replacement of sums by integrals in the infinite volume limit, we find 



£:(0)= lim^ = i02(M2-A02 



/.^<A^(^v/^^+M^-!Ab 



4>^ + Id{M\A) 



where M = M[cj)) depends on (j) through the gap equation ([4.17|) . The explicit calculation 
of the integrals involved shows that the energy density difference ^(0) — ^(0) [which for 
unbroken symmetry is nothing but the effective potential V^ff(0)], is indeed finite in the 
limit A — > oo, as required by a correct renormalization scheme. Notice that the finiteness of 
the energy density difference can be showed also by a simpler and more elegant argument, 
as presented below in section |1V IJ| . This check would fail instead when D = 3 if only the 
bare coupling constant Ab would appear in the last formula. 

The tdHF approximation derived above represents a huge simplification with respect to 
the original problem, but its exact solution still poses itself as a considerable challenge. As a 
matter of fact, a numerical approach is perfectly possible within the capabilities of modern 
computers, provided the number of equations ( [4.5| ) is kept in the range of few thousands. 
As will become clear later on, even this numerical workout will turn out not to be really 
necessary in the form just alluded to, at least for the purposes of this paper. 



A. On symmetry breaking 

Quite obviously, in a finite volume and with a UV cutoff there cannot be any symmetry 
breaking, since the ground state is necessarily unique and symmetric when the number of 
degrees of freedom is finite |^. However, we may handily envisage the situation which 



would imply symmetry breaking when the volume diverges. 
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Let us first consider the case that we would call of unbroken symmetry. In this case 
the HF ground state is very close to the member with = of the family of gaussian 
states introduced before. The difference is entirely due to the quartic terms in H},. This 
correction vanish when L -^ oo, since all wavefunctions ipk have L— independent widths, 
so that one directly obtains the symmetric vacuum state with all the right properties of 
the vacuum (translation invariance, unicity, etc.) upon which a standard scalar massive 
particle Fock space can be based. The HF approximation then turns out to be equivalent 
to the resummation of all "cactus diagrams" for the particle self-energy |p. In a finite 
volume, the crucial property of this symmetric vacuum is that all frequencies uj\ are strictly 
positive. The generalization to non-equilibrium initial states with 7^ is rather trivial: 
it amounts to a shift by L^/'^cj) on V^olv^o)- In the limit L — i> 00 we should express iI)q as 
a function of ^ = L'^^'^^po so that, iV'olOP ~^ ^(^ ~ 0); while all other wavefunctions ipk 
will reconstruct the gaussian wavefunctional corresponding to the vacuum |0, M) of a free 
massive scalar theory whose mass M = M{(f)) solves the gap equation ( 4.17| ). The absence 



oi ipQ in |0, M) is irrelevant in the infinite volume limit, since ((/Jq) = L^(f)'^+ terms of order 
L°. The effective potential Ves{4>) = S{(i>) — S{0), where £{4>) is the lowest energy density 
at fixed (f) and infinite volume, is manifestly a convex function with a unique minimum in 
= 0. 

Now let us consider a different situation in which one or more of the ul are negative. 
Quite evidently, this might happen only for k small enough, due to the /c^ in the gap equation 
[thus eq. ([4.2CI|) remains valid and the ultraviolet renormalization is the same as for unbroken 
symmetry]. Actually we assume here that only cJq < 0, postponing the general analysis. 
Now the quartic term in Hq cannot be neglected as L — > 00, since in the ground state ipo is 
symmetrically concentrated around the two minima of the potential ^uJq^Pq + -^^^Pq, that is 
ipo = ±{~uJqL^ /Xy/"^. If we scale ipo as ipo = L^^'^i then Hq becomes 

H. = ~^t+'^Ue + -A (4.21) 



so that the larger L grows the narrower ipo{C) becomes around the two minima ^ = 
±(— cJo/A)^/^. In particular {S,"^) -^ —^l/X when L —* 00 and (v^q) ^ L^{^'^). More- 
over, the energy gap between the ground state of Hq and its first, odd excited state as well 
as difference between the relative probability distributions for ^ vanish exponentially fast in 
the volume L^. 

Since by hypothesis all ul with k ^ are strictly positive, the ground state 'ipk with 
k ^ are asymptotically gaussian when L ^ 00 and the relations ( [4.11|) tend to the form 

u;l = e + M^ = e + m^ 
M^ = -2ujI = ml + 3K{L-^{ipl) + Sq) = ml + SX^^uj^q + SAb/^lm^, A)] 

This implies the identification Uq = —m^jl and the bare mass parameterization 

ml=(\- |Ab/A) m" - 'iX^,Ijj{jr? , A) (4.22) 



characteristic of a negative cJq [compare to eq. ( |4.16| )], with m the physical equilibrium 



mass of the scalar particle, as in the unbroken symmetry case. The coupling constant 
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renorinalization is the same as in eq. ( |4.18| ) as may be verified by generalizing to the 



minimum energy states with given field expectation value 0; this minimum energy is nothing 
but the effective potential Veflf(0); of course, since ipo is no longer asymptotically gaussian, 
we cannot simply shift it by L^^'^cf) but, due to the concentration of iI)q on classical minima 
as L ^ oo, one readily finds that 14^(0) is the convex envelope of the classical potential, 
that is its Maxwell construction. Hence we find 




-^0 



and the gap equation for the 0— dependent mass M can be written, in terms of the step 
function B and the extremal ground state field expectation value v = r7i/v2A, 

M^ = m'^ + 3Ab(0^ - v^) 6(02 - v^) + 3Ab [/d(M^ A) - /^(m^ A)] (4.23) 



We see that the specific bare mass parameterization ( [4.22|) guarantees the non- 



renormalization of the tree-level relation f ^ = vr? j2\ ensuing from the typical symmetry 
breaking classical potential V[<\)) = ^\{4>'^—v^)'^. With the same finite part prescription as in 
eq. (|4.171) , the gap equation ([4.23|) leads to the standard coupling constant renormalization 



(|A|) when D = 3. 

In terms of the probability distributions |V'o(OP ^^ the scaled amplitude ^ = L~^/'^(fo, 
the Maxwell construction corresponds to the limiting form 

L^oo I 5(^ - 0) , 0^ > -U^ 

On the other hand, if Uq is indeed the only negative squared frequency, the k j^ part of 
this minimum energy state with arbitrary = {4>{x)) is better and better approximated 
as A — >■ cxD by the same gaussian state |0, M) of the unbroken symmetry state. Only the 
effective mass M has a different dependence M{(j)), as given by the gap equation ( [4.23[ ) 
proper of broken symmetry. 

As in the large A^ approach, at infinite volume we may write 



2VWTW 



where C(0) = 0^ in case of unbroken symmetry (that is Wg > 0), while C(0) = max(f 2, 0^) 
,2 




when cUq < 0. This corresponds to the field correlation in space 



(0(x)0(y)) = / ^(^De*(^-^) = C7(0) + Az,(x - y, M) 

where ^d{x ~ UiM) is the massive free field equal-time two points function in D space 
dimensions, with self-consistent mass M. As in large A^, the requirement of clustering 



(0(x)0(y)) ^ (0(x))2 = V 
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contradicts the infinite volume limit of 

k 

except at the (now only two) extremal points = ±v. In fact we know that the L ^ oo 
limit of the finite volume states with cfp' < v"^ violate clustering, because the two peaks of 
iPq{C} have vanishing overlap in the limit and the first excited state becomes degenerate 
with the vacuum: this implies that the relative Hilbert space splits into two orthogonal Fock 
sectors each exhibiting symmetry breaking, (0(x)) = ±f, and corresponding to the two 
independent equal weight linear combinations of the two degenerate vacuum states. The 
true vacuum is either one of these symmetry broken states. Since the two Fock sectors are 
not only orthogonal, but also superselected (no local observable interpolates between them), 
linear combinations of any pair of vectors from the two sectors are not distinguishable from 
mixtures of states and clustering cannot hold in non-pure phases. It is perhaps worth 
noticing also that the Maxwell construction for the effective potential, in the infinite volume 
limit, is just a straightforward manifestation of this fact and holds true, as such, beyond the 
HF approximation. 

To further clarify this point and in view of subsequent applications, let us consider the 
probability distribution for the smeared field 0/ = J d^x (l){x)f{x), where 



/(-) = /(-•) = ^ E A- -■'■' ,~„ / 75^ m e-- 



is a smooth real function with / d^x f{x) = 1 {i.e. /o = 1) localized around the origin (which 
is good as any other point owing to translation invariance). Neglecting in the infinite volume 
limit the quartic corrections for all modes with A; 7^ 0, so that the corresponding ground 
state wavefunctions are asymptotically gaussian, this probability distribution evaluates to 



(-^y /•+00 

Vi{u<(f)f<u + du) = ---—-j:r d^|^o(OI exp 

ZTTZjf r'^ J-00 



where 






d^k f{kf 



(27r)^ 2v/FT^ 



In the unbroken symmetry case we have |'0o(OP ~ ^{^ — 0) as L — i> 00, while the limiting 
form ( |4.24| ) holds for broken symmetry. Thus we obtain 



Fr{u<(f)f<u + du)=pf{u — (f))du, P/(m) = (27rSj) ' exp 



2' 



2S, 



for unbroken symmetry and 



\pf{u — (p) du , 0^ > f^ 
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for broken symmetry. Notice that the momentum integration in the expression for Sj needs 
no longer an ultraviolet cutoff; of course in the limit of delta-like test function f{x), Hf 
diverges and Pf{u) flattens down to zero. The important observation is that Pr(M < (pf < 
u + du) has always a single peak centered in m = for unbroken symmetry, while for broken 
symmetry it shows two peaks for 0^ < v"^ and Sj small enough. For instance, if = 0, then 
there are two peaks for Sj < v"^ [implying that f{k) has a significant support only up to 
wavevector k of order f , when D = 3, or r77,exp(const v"^) when D = 1]. 

To end the discussion on symmetry breaking, we may now verify the validity of the 
assumption that only uJq is negative. In fact, to any squared frequency ul (with k ^ 0) that 
stays strictly negative as L ^ oo there corresponds a wavefunction ipk that concentrates on 
(fl + ip'^_f^ = —ulL^/X ; then eqs. (|4.11| ) implies —2ul = k^ + rn? for such frequencies, while 
ujI = k"^ + m? for all frequencies with positive squares; if there is a macroscopic number of 
negative ujI (that is a number of order L^)., then the expression for uJq in eq. ( [4.11|) will 



contain a positive term of order L^ in the r.h.s., clearly incompatible with the requirements 
that cJq < and m^ be independent of L; if the number of negative ul is not macroscopic, 
then the largest wavevector with a negative squared frequency tends to zero as L ^ oo (the 
negative ujI clearly pile in the infrared) and the situation is equivalent, if not identical, to 
that discussed above with only oJq < 0. 

B. Out of equilibrium dynamics 

We considered above the lowest energy states with a predefinite uniform field expectation 
value, {4>{x)) = 0, and established how they drastically simplify in the infinite volume limit. 
For generic these states are not stationary and will evolve in time. By hypothesis V'fc is 
the ground state eigenfunction of H^ when k > 0, and therefore \ipk\'^ would be stationary 
for constant Uk, but ipo is not an eigenfunction of Hq unless = 0. As soon as \ipk\'^ starts 
changing, {ipD changes and so do all frequencies Uk which are coupled to it by the eqs. 
( |4.11| ). Thus the change propagates to all wavefunctions. The difficult task of studying 
this dynamics can be simplified with the following scheme, that we might call gaussian 
approximation. We first describe it and discuss its validity later on. 

Let us assume the usual gaussian form for the initial state [see eq. ( |4.13| ) and the 
discussion following it]. We know that it is a good approximation to the lowest energy state 
with given (ipo) for unbroken symmetry, while it fails to be so for broken symmetry, only as 
far as ipo is concerned, unless 0^ > v"^. At any rate this is an acceptable initial state: the 
question is about its time evolution. Suppose we adopt the harmonic approximation for all 
Hk with /c > by dropping the quartic term. This approximation will turn out to be valid 
only if the width of ipk do not grow up to the order L^ (by symmetry the center will stay 
in the origin). In practice we are now dealing with a collection of harmonic oscillators with 
time-dependent frequencies and the treatment is quite elementary: consider the simplest 
example of one quantum degree of freedom described by the gaussian wavefunction 



^^^''^=(2;^^^'^P 



1/1 .S^ 2 



2 V2(t2 a 



i-\q 



where s and a are time-dependent. If the dynamics is determined by the time-dependent 
harmonic hamiltonian ^[—5^ + ^^(t)^^^], then the Schroedinger equation is solved exactly 
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provided that s and a satisfy the classical Hamilton equations 

1 



a = s , s = —up' a + 



4(t3 
It is not difficult to trace the "centrifugal" force (4o")^'^ which prevent the vanishing of a to 



Heisenberg uncertainty principle [Il2| , p!6 



The extension to our case with many degrees of freedom is straightforward and we find 
the following system of equations 

2-7/^0 = a^^^ , -^ = -u;lak + ^, k>0 (4.25) 



coupled in a mean-field way by the relations ( f4.11| ), which now read 



ul = k^ + M^- 6A L-^al , k>0 



M' = ml + 3Ab {l-^{vI) + So) , So = -^ ^ 



1 ^2 (4.26) 

fD Z^^k 



This stage of a truly quantum zero-mode and classical modes with A; > does not appear 
fully consistent, since for large volumes some type of classical or gaussian approximation 
should be considered for ipo too. We may proceed in two (soon to be proven equivalent) 

ways: 

1. We shift ifo = L^/'^cf) + 770 and then deal with the quantum mode ?7o in the gaussian 
approximation, taking into account that we must have (r/o) = at all times. This is 
most easily accomplished in the Heisenberg picture rather than in the Schroedinger 
one adopted above. In any case we find that the quantum dynamics of y^o is equivalent 
to the classical dynamics of and o"o = iVo)^^'^ described by the ordinary differential 
equations 



— = -u;,c^-X^ , ^ = --0-0 + ^3 



-^o<^-A0^ -irf = -^0 ^0 + 7-3 (4.27) 







where u^ = M' - 3A L-^((/?^) and (y?^) = L^0^ + a, 



2. We rescale (^0 = L^^"^^ right away, so that Hq takes the form of eq. ([4.21|) . Then L 



cxD is the classical limit such that '?/'o(0 concentrates on ^ = which evolves according 
to the ffist of the classical equations in ( [4.27| ). Since now there is no width associated 



to the zero-mode, is coupled only to the widths ak with k ^ hj Uq = M^ — 3A0^, 
while M'^ = ml + 3Ab(02 + So). 

It is quite evident that these two approaches are completely equivalent in the infinite 
volume limit, and both are good approximation to the original tdHF Schroedinger equations, 
at least provided that a^ stays such that L^^ctq vanishes in the limit for any time. In this 
case we have the evolution equations 



1^.PA^=-M^)^, if^.-(,= ^MV..jl3 



(2A0^-M^)0, -^ = -(k' + M')ak + ^r^ (4.28) 

k 
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mean-field coupled by the L ^ oo limit of eqs. ( [4.26| ), namely 



M' 



w? + 3Ab 



T.-Io{m\A) 



for unbroken symmetry [that is rn^ as in eq. ( [4.16| )] or 

M^ = m^ + 3Ab [02 - t;2 + S - /^(m^ A) 



m 



2\v' 



(4.29) 



(4.30) 



for broken symmetry [that is m^ as in eq. ([4.22|) ]. In any case we define 



TdE^, 



LD 



d^k 



oo Jk'2<A'2 (27r) 



D 



CTi^ 



as the sum, or integral, over all microscopic gaussian widths [N.B.:this definition differs from 
that given before in eq. ( [4.8|) by the classical term 0^]. 

The conserved HF energy (density) corresponding to these equations of motion reads 



V 



^^^)'+2l^ 



E^ 



CTl, 



Mi' 



(4.31) 



2 ^2 I 



W, 



+ im2(02 + s) + |Ab(02 + S)2 - iA0^ 



Up to additive constants and terms vanishing in the infinite volume limit, this expression 
agrees with the general HF energy of eq. ([4.12| ) for gaussian wavef unctions. It holds both 
for unbroken and broken symmetry, the only difference being in the parameterization of the 
bare mass in terms of UV cutoff and physical mass, eqs. ( [4.16| ) and ( [4.22| ). The similarity 
to the energy functional of the large N approach, eq. ( [4.31|) , is evident; the only difference, 
apart from the obvious fact that ^ is a 0{n) vector rather than a single scalar, is in the 
mean-field coupling ak-fj) and Ufc-S, due to different coupling strength of transverse and 
longitudinal modes. 

This difference between the HF approach for discrete symmetry {i.e N = 1) and the 
large A^ method for the continuous 0(A^)-symmetry is not very relevant if the symmetry is 
unbroken [it does imply however a significantly slower dissipation to the modes of the back- 
ground energy density]. On the other hand it has a drastic consequence on the equilibrium 
properties and on the out-of-equilibrium dynamics in case of broken symmetry (see below), 
since massless Goldstone bosons appear in the large N approach, while the HF treatment 
of the discrete symmetry case must exhibits a mass also in the broken symmetry phase. 

The analysis of physically viable initial conditions proceeds exactly as in the large N 
approach and will not be repeated here, except for an important observation in case of 
broken symmetry. The formal energy minimization w.r.t. ak at fixed leads again to eqs. 
( ^.18| ), and again these are acceptable initial conditions only if the gap equation that follows 
from eq. ([4.30|) in the L ^ oo limit, namely 



M' 



m^ + 3Ab 02 - t;2 + Id{M\ A) - /^(m^. A) 



(4.32) 



admits a nonnegative, physical solution for M^. Notice that there is no step function in eq. 
( [4.32| ), unlike the static case of eq. (|4.23|) , because Uq was assumed to be microscopic, so 
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that the infinite volume a1 has no delta-hke singularity in A; = 0. Hence M = m solves 
eq. (^4.32 ) only at the extremal points (j) = ±f , while it was the solution of the static gap 



equation ( [4.23| ) throughout the Maxwell region —v < (j) < v. The important observation is 
that eq. ([4.32|) admits a positive solution for M^ also within the Maxwell region. In fact it 
can be written, neglecting as usual the inverse-power corrections in the UV cutoff 



A(M) A 

and there exists indeed a positive solution M^ smoothly connected to the ground state, 
0^ = f^ and M^ = rn^, whenever 0^ > v'^/3. The two intervals v"^ > cjP' > f^/3 correspond 
indeed to the metastability regions, while (p < v'^/3 is the spinodal region, associated to a 
classical potential proportional to (0^ — f ^)^. This is another effect of the different coupling 
of transverse and longitudinal modes: in the large N approach there are no metastability 
regions and the spinodal region coincides with the Maxwell one. As in the large A^ approach 
in the spinodal interval there is no energy minimization possible, at fixed background and 
for microscopic widths, so that a modified form of the gap equation like eq. (|3.2CI|) should 
be applied to determine ultraviolet-finite initial conditions. 

The main question now is: how will the gaussian widths 0"^ grow with time, and in 
particular how will (Tq grow in case of method 1 above, when we start from initial conditions 
where all widths are microscopic? For the gaussian approximation to remain valid through 
time, all ak, and in particular ctq, must at least not become macroscopic. In fact we have 
already positively answered this question in the large N approach and the HF equations 
( [4.28| ) do not differ so much to expect the contrary now. In particular, if we consider the 

special initial condition = = 0, the dynamics of the widths is identical to that in the 
large A^ approach, apart from the rescaling by a factor of three of the coupling constant. 
Thus our HF approximation confirms the large A^ approach in the following sense: even if 
one considers in the variational ansatz the possibility of non-gaussian wavefunctionals, the 
time evolution from gaussian and microscopic initial conditions is effectively restricted for 
large volumes to non-macroscopic gaussians. 

Strictly speaking, however, this might well not be enough, since the infrared fluctuations 
do grow beyond the microscopic size to become of order L. Then the quartic term in the 
low— fc Hamiltonians Hk is of order L and therefore it is not negligible by itself in the 
L ^ oo limit, but only when compared to the quadratic term, which for a fixed ujI of order 
1 would be of order L^. But we know that, when = 0, after the spinodal time and before 
the Ti, the effective squared mass M^ oscillates around zero with amplitude decreasing as 
t~^ and a frequency flxed by the largest spinodal wavevector. In practice it is "zero on 
average" and this reflect itself in the average linear growth of the zero-mode fluctuations 
and, more generally, in the average harmonic motion of the other widths with non-zero 
wavevectors. In particular the modes with small wavevectors of order L~^ feel an average 
harmonic potential with a;^ of order L~^ . This completely compensate the amplitude of the 
mode itself, so that the quadratic term in the low— /c Hamiltonians Hk is of order L", much 
smaller than the quartic term that was neglected beforehand in the gaussians approximation. 
Clearly the approximation itself no longer appear fully justifled and a more delicate analysis 
is required. We intend to return on this issue in a future work, restricting in the next section 
our observations on some rather peculiar consequences of the gaussians approximation that 
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provide further evidence for its internal inconsistency. 



C. Effective potential and late time evolution 

By definition, the gaussian approximation of the effective potential Kflf(0) coincides with 
the infinite-volume limit of the potential energy V(0, {(Tk}) of eq. ( 4.31J ) when the widths are 



of the ^—dependent, energy-minimizing form ( |3.18| ) with the gap equation for M^ admitting 
a nonnegative solution. As we have seen, this holds true in the unbroken symmetry case for 
any value of the background 0, so that the gaussian V^g is identical to the HF one, since 
all wavefunctions ipk are asymptotically gaussians as L — >■ oo. In the presence of symmetry 
breaking instead, this agreement holds true only for 0^ > f^; for f^/3 < 0^ < f^ the 
gaussian V^s exists but is larger than the HF potential, which is already fiat. In fact, for 
any 0^ > f^/3, we may write the gaussian Ves as 

V,si^) = Kff(-0) = V,siv) + /"^ duu[M{uf - 2Am2] 

Jv 

where M(m)^ solves the gap equation (|4.33|) , namely M(m)^ = X{M{u)){3u'^ — v"^). In each 
of the two disjoint regions of definition this potential is smooth and convex, with unique 
minima in +v and — f , respectively. Its HF counterpart is identical for cf? > v"^, while it 
takes the constant value Ves{v) throughout the internal region 0^ < f^. On the other hand 
the gaussian Ves cannot be defined in the spinodal region 0^ < v^/3, where the gap equation 
does not admit a nonnegative solution in the physical region far away from the Landau pole. 
Let us first compare this HF situation with that of large A^. There the different coupling 
of the transverse modes, three time smaller than the HF longitudinal coupling, has two main 
consequences at the static level: the gap equation ( |3l8| ) does not admit nonnegative solutions 
for (j) < f^, so that the spinodal region coincides with the region in which the effective 
potential is flat, and the physical mass vanishes. The out-of-equilibrium counterpart of 
this is the dynamical Maxwell construction: when the initial conditions are such that (p 
has a limit for t — > oo, all possible asymptotic values exactly span the flatness region (and 
the effective mass vanishes in the limit). In practice this means that |0| is not the true 
dynamical order parameter, whose large time limit coincides with v, the equilibrium fleld 
expectation value in a pure phase. Rather, one should consider as order parameter the 
renormalized local (squared) width 

where the last equality follows from the deflnition itself of the effective mass M [see eq. 



-,2 



( p.l7| )]. Since M vanishes as t ^ oo when tends to a limit within the flatness region. 



we flnd the renormalized local width tends to the correct value v which characterizes the 
broken symmetry phase, that is the bottom of the classical potential. 

Let us now examine what happens instead in the dynamics of the HF approximation, 
where at the static level the spinodal region 0^ < v^/3 is smaller than the flatness region 
(fp' < v"^. Our (preliminary) numerical evidence shows that a dynamical Maxwell construction 
take place as in the large N approximation, but it covers only the spinodal region [see for 
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instance figs. |^ and |T0|. If tlie background starts witli zero velocity inside tlie spinodal 
interval, then it tends to a limit within the same interval, the asymptotic force vanishing 
because M^ = 2A0^ in the limit [see eq. ( [4.28| )]. Hence we find that the order parameter 



/./ ^2^ 72 ^ ^' M2-2A, 



3 3A 

[the relation (|4.30|) was used in the last equality] tends to v'^/3 not v"^. It "stops at the 
spinodal line". This fact is at the basis of the so-called spinodal inflation [^]. Even without 
a conspicuous numerical evidence for the full dynamical Maxwell construction, this results is 
manifestly true when initially (and therefore at any time) = = 0, because the equations 
of motion for the widths are the same as in large N and the effective mass only differs 
by the factor of three in front of the quantum back-reaction S. Since practically all the 
back-reaction takes place during the exponential growth of the unstable spinodal modes, 
one could say that, in the gaussian HF approximation, this back-reaction is to strong and 
prevents the quantum fluctuations from sampling the true minima of the classical potential. 
A side effect of this is that the effective mass M does not tend to its correct equilibrium 
value m as t — > oo, unlike in the large N approach. 

In the previous section we have discussed the possible origin of the problem: in our tdHF 
approach the initial gaussian wavefunctions are allowed to evolve into non-gaussian forms, 
but they simply do not do it in a macroscopic way, within a further harmonic approxima- 
tion for the evolution, so that in the infinite-volume limit they are indistinguishable from 
gaussians at all times. But when M^ is on average not or order L°, but much less, as it 
happens for suitable initial conditions, infrared modes of order L will be dominated by the 
quartic term in our Schroedinger equations ( [4.5D , showing a possible internal inconsistency 
of the gaussians approximation. 

Another interesting point concerns the dynamical Maxwell construction itself, within the 
gaussians approximation. In fact it is not at all trivial that the effective potential, in any of 
the approximation previously discussed, does bear relevance on the asymptotic behavior of 
the infinite-volume system whenever a fixed point is approached. Strictly speaking in fact, 
even in such a special case the effective potential has little to say on the dynamics, since it is 
obtained from a static minimization and the energy is not at its minimum at the initial time 
and is exactly conserved in the evolution. On the other hand, if a solution of the equations of 
motion ( [4.28| ) exists in which the background tends to a constant 0oo as t ^ cxd, one might 
expect that the effective action (which however is nonlocal in time) somehow reduces to a 
(infinite) multiple the effective potential, so that 0oo should be an extremal of the effective 
potential. This is still an open question that deserve further analytic studies and perhaps 
some further numerical confirmation. 

V. CONCLUSIONS AND PERSPECTIVES 

In this work we have presented a rather detailed study of the dynamical evolution out of 
equilibrium, in finite volume (a cubic box of size L in 3D), for the 0^ QFT. For comparison, 
we have also analyzed some static characteristics of the theory both in unbroken and broken 
symmetry phases. We have worked in two non-perturbative approximation schemes, namely 
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the large A^ expansion in the leading order and a generalized time dependent Hartree-Fock 
approximation. 

We have reached two main conclusions. 

The first, based on strong numerical evidence, is that the linear growth of the zero-mode 
quantum fluctuations, observed already in the large N approach of refs. P-p!0|], cannot be 
consistently interpreted as a new form of Bose-Einstein condensation. In fact, in finite 
volume, (To never grows to 0(L^/^) if it starts from a microscopic value, that is at most of 
order L^/^. On the other hand all long-wavelength fluctuations rapidly become of order L, 
signalling a novel infrared behavior quite different from free massless fields at equilibrium 
[recall that the large A^ or HF approximations are of mean field type, with no direct interac- 
tion among particle excitations]. This is in agreement with the properties of the two-point 
function determined in [Q. 

The second point concerns the HF approach to the out-of-equilibrium 0^ QFT. We 
have shown that, within a slightly enlarged tdHF approach that allows for non-gaussian 
wavefunctions, one might recover the usual gaussians HF approximation in a more controlled 
way, realizing that the growth of long-wavelength fluctuations to order L in fact undermines 
the self-consistency of the gaussians approximation itself. The first manifestation of this 
weakness is the curious "stopping at the spinodal line" of the width of the gaussian quantum 
fluctuations. This does not happen in the large N approach because of different coupling 
of transverse mode (the only ones that survive in the N ^ oo limit) with respect to the 
longitudinal modes of the A^ = 1 case in the HF approach. 

Clearly further study, both analytical and numerical, is needed in our tdHF approach to 
better understand the dynamical evolution of quantum fluctuations in the broken symmetry 
phase. Another interesting direction is the investigation, in our HF approximation, of the 
case of finite A^, in order to interpolate smoothly the results for A^ = 1 to those of the 1/A^ 
approach. One first question concerns whether or not the Goldstone theorem is respected in 
the HF approximation [^. It would be interesting also to study the dynamical realization 



of the Goldstone paradigm, namely the asymptotic vanishing of the effective mass in the 
broken symmetry phases, in different models; this issue needs further study in the 2D case 



r^ , where it is known that the Goldstone theorem is not valid. 

Another open question concerns the connection between the minima of the effective 
potential and the asymptotic values for the evolution of the background, within the simplest 
gaussian approximation. As already pointed out in |^, a dynamical Maxwell construction 
occurs for the 0{N) model in infinite volume and at leading order in 1/A^ in case of broken 
symmetry, in the sense that any value of the background within the spinodal region can 
be obtained as large time limit of the evolution starting from suitable initial conditions. 
Preliminary numerical evidence suggests that something similar occurs also in the Hartree 
approximation for a single field, but a more thorough and detailed analysis is needed. 
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APPENDIX A: DETAILS OF THE NUMERICAL ANALYSIS 



We present here the precise form of the evolution equations for the field background and 
the quantum mode widths, which control the out-of-equilibrium dynamics of the 0^ model 
in finite volume at the leading order in the 1/A^ approach and in the tdHF approach, as 
described in sections [III CI and rVBl. We restrict here our attention to the tridimensional 



case. 



Let us begin by noticing that each eigenvalue of the Laplacian operator in a 3D finite 
volume is of the form ^n = (t") '^' "^here n is a non-negative integer obtained as the sum 
of three squared integers, n 
each eigenvalue, representing the number of different ordered triples {n^ 



nl + nl + nl- Then we associate a degeneracy factor (?„ to 



"j/; "^z 



yielding the 
same n. One may verify that Qn takes on the continuum value of Aiik'^ in the infinite volume 

limit, where k = (^j n is kept fixed when L -^ oo. 

Now, the system of coupled ordinary differential equations is, in case of the large N 
approach. 



•P .,2 










0"n — 



4a3 







(Al) 



where the index n ranges from to A/"^, A/" = AL/27r and M'^{t) is defined by the eq. 
( p.l6| ) in case of unbroken symmetry and by eq. ( p. 171) in case of broken symmetry. The 
back-reaction S reads, in the notations of this appendix 



A^2 



L^ 



T.9n 
n=0 



cr„ 



Technically it is simpler to treat an equivalent set of equations, which are formally linear 
and do not contain the singular Heisenberg term oc a~^. This is done by introducing the 
complex mode amplitudes z^ = 0"„exp(i6'„), where the phases 9n satisfy <j'^9n = 1. Then we 
find a discrete version of the equations studied for instance in ref fl^ , 



d"" 



+ 



27r 



subject to the Wronskian condition 



n + M^ 



N^ 







L^ 



/ , 9n\Zn\ 



(A2) 



n=0 



Zn, Zn 



Zn. Zn 



One realizes that the Heisenberg term in cr„ corresponds to the centrifugal potential for the 
motion in the complex plane of z„. 

Let us now come back to the equations (^). To solve these evolution equations, we 
have to choose suitable initial conditions respecting the Wronskian condition. In case of 
unbroken symmetry, once we have fixed the value of <^ and its first time derivative at initial 
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time, the most natural way of fixing the initial conditions for the Zn is to require that they 
minimize the energy at t = 0. We can obviously fix the arbitrary phase in such a way to 
have a real initial value for the complex mode functions 

2„(0) 



V2fin 




where Vtn = Jk^ + M^(0). The initial squared effective mass M^(t = 0), has to be deter- 
mined self-consistently, by means of its definition ( |3.16| ). 



In case of broken symmetry, the gap equation is a viable mean for fixing the initial 
conditions only when lies outside the spinodal region [cfr. eq (|3.21|)]; otherwise, the gap 



equation does not admit a positive solution for the squared effective mass. In that case, 
we have to resort to other methods, in order to choose the initial conditions. Following the 
discussion presented in [IIIC] , one possible choice is to set al = — ^= — = for k^ < \M^\ in 



eq. ( |3.2CI|) and then solve the corresponding gap equation ( p.2(]| ). An other acceptable choice 



would be to solve the gap equation (|3.2CI|) , once we have set a massless spectrum for all the 
spinodal modes but the zero mode, which is started from an arbitrary, albeit microscopic, 
value. 

There is actually a third possibility, that is in some sense half a way between the unbroken 
and broken symmetry case. We could allow for a time dependent bare mass, in such a way 
to simulate a sort of cooling down of the system. In order to do that, we could start with a 
unbroken symmetry bare potential (which fixes initial conditions naturally via the gap equa- 
tion) and then turn to a broken symmetry one after a short interval of time. This evolution 
is achieved by a proper interpolation in time of the two inequivalent parameterizations of 
the bare mass, eqs. (|0| ) and (|3.10|) . 

We looked for the influence this different choices could produce in the results and indeed 
they depend very little and only quantitatively from the choice of initial condition we make. 

As far as the numerical algorithm is concerned, we used a 4th order Runge-Kutta algo- 
rithm to solve the coupled differential equations (|A^), performing the computations in boxes 
of linear size ranging from L = 10 to L = 200 and verifying the conservation of the Wron- 
skian to order 10~^. Typically, we have chosen values of A/" corresponding to the UV cutoff 
A equal to small multiples of m for unbroken symmetry and of vyX for broken symmetry. 
In fact, the dynamics is very weakly sensitive to the presence of the ultraviolet modes, once 
the proper subtractions are performed. This is because only the modes inside the unstable 
(forbidden or spinodal) band grow exponentially fast, reaching soon non perturbative ampli- 
tudes (i.e. ~ A"^/^), while the modes lying outside the unstable band remains perturbative, 
contributing very little to the quantum back-reaction and weakly affecting the overall 
dynamics. The unique precaution to take is that the initial conditions be such that the 
unstable band lay well within the cutoff. 
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FIGURES 
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FIG. 1. Zero-mode amplitude evolution for different values of the size L = 20,40,60,80,100, 
for A = 0.1 and broken symmetry, with (f) = 0. 




FIG. 2. Time evolution of the squared effective mass M^ in broken symmetry, for L = 100 and 
A = 0.1. 
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FIG. 3. The quantum hack-reaction S, with the parameters as in Fig. 
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FIG. 4. Zero-mode amplitude evolution for different value of the renormalized coupling constant 
A = 0.01, 0.1, 1, for L = 100 and broken symmetry, with (j) = 0. 
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FIG. 5. Detail of M^ near t = t^ for L = 40 (dotted line). The case L = 80 is plotted for 
comparison (solid line). 
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FIG. 6. Detail of ctq near t = tl for L = 40 (dotted line). The case L = 80 is plotted for 
comparison (solid line). 
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FIG. 7. Next-to-zero mode (k = 271 /L) amplitude evolution for different values of the size 
L = 20,40, 60, 80, 100, for A = 0.1 and broken symmetry, with (/) = 0. 




100 200 300 400 500 600 700 800 900 1000 

FIG. 8. Next-to-zero mode (k = 2tt/L) amplitude evolution for different value of the renor- 
malized coupling constant A = 0.01, 0.1, 1, for L = 100 and broken symmetry, with (p = 0. 
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FIG. 9. Evolution of the background for two different initial conditions within the spinodal 
interval, in the tdHF approximation: (j){t = 0) = 0.1 (dotted line) and (j)(t = 0) = 0.4 (solid line). 




FIG. 10. Evolution of M^ for the two initial conditions of fig. |^. 
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